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ABSTRACT 

We present near-infrared coronagraphic imaging polarimetry of RY Tau. The scattered light 
in the circumstellar environment was imaged at i7-band at a high resolution (~0”.05) for the 
first time, using Subaru-HiCIAO. The observed polarized intensity (PI) distribution shows a 
butterfly- like distribution of bright emission with an angular scale similar to the disk observed 
at millimeter wavelengths. This distribution is offset toward the blueshifted jet, indicating the 
presence of a geometrically thick disk or a remnant envelope, and therefore the earliest stage of 
the Class II evolutionary phase. We perform comparisons between the observed PI distribution 
and disk models with: (1) full radiative transfer code, using the spectral energy distribution 
(SED) to constrain the disk parameters; and (2) monochromatic simulations of scattered light 
which explore a wide range of parameters space to constrain the disk and dust parameters. We 
show that these models cannot consistently explain the observed PI distribution, SED, and the 
viewing angle inferred by millimeter interferometry. We suggest that the scattered light in the 
near-infrared is associated with an optically thin and geometrically thick layer above the disk 
surface, with the surface responsible for the infrared SED. Half of the scattered light and thermal 
radiation in this layer illuminates the disk surface, and this process may significantly affect the 
thermal structure of the disk. 

Subject headings : protoplanetary disks — stars: indivdual (RY Tau) — stars: pre-main sequence — 
polarization 


1 Institute of Astronomy and Astrophysics, Academia 
Sinica. P.O. Box 23-141, Taipei 10617, Taiwan, R.O.C.; 
hiro@asiaa.sinica.edu. tw 

2 National Astronomical Observatory of Japan, 2-21-1 
Osawa, Mitaka, Tokyo 181-8588, Japan 

3 H.L. Dodge Department of Physics and Astronomy, 


University of Oklahoma, 440 W Brooks St Norman, OK 

73019, USA 

4 Max Planck Institute for Astronomy, Koenigstuhl 17, 
D-69117 Heidelberg, Germany 

5 Eureka Scientific, 2452 Delmer, Suite 100, Oakland CA 
96002, USA 

6 Institute for Astronomy, University of Hawaii, 640 


1 


1. Introduction 

Optical and near-infrared observations have re- 
vealed structures in protoplanetary disks at the 
highest angular resolutions currently available. 
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This has provided powerful tools for investigat- 
ing the possibility of ongoing planet formation 
and to test the related theories. In particular, 
the technique of coronagraphic imaging has been 
extensively used to suppress the stellar flux and 
detect scattered light from the disk surface with 
high sensitivities (see Watson et al. 2007, for a 
review). Without this technique, resolved images 
of protoplanetary disks can only be obtained in 
limited circumstences at optical and near-infrared 
wavelengths, via an edge-on view or silhouette 
against bright background nebular emission (e.g., 
McCaughrean et al. 2000, for a review). 

Using coronagraphy, the ongoing survey “Strate- 
gic Explorations of Exoplanets and Disks with 
Subaru” (SEEDS, Tamura 2009) with Subaru- 
HiCIAO (Tamura et al. 2006) and AO 188 (Hayano 
et al. 2004) has recently discovered structures in a 
number of protoplanetary disks (Thalmann et al. 
2010), in particular in the polarized intensity 
(hereafter PI) distribution in the near infrared 
(Hashimoto et al. 2011, 2012; Muto et al. 2012; 
Kusakabe et al. 2012; Tanii et al. 2012; Mayama 
et al. 2012; Dong et al. 2012a; Grady et al. 2013; 
Follette et al. 2013). PI imaging has been used 
for observations of most of the disks in the SEEDS 
program, since this suffers significantly less con- 
tamination from the stellar flux than the normal 
intensity /. Some disks are associated with spi- 
ral structures, disk holes, or azimuthal gaps in 
ring-like flux distributions, which are potential 
signatures of ongoing planet formation. The scat- 
tered light from the disk is also useful for probing 
grain growth (e.g., Min et al. 2012; Tanii et al. 
2012) which may be related to formation of rocky 
cores. 

The goal of the SEEDS program for protoplan- 
etary disks is to observe a large number of objects 
with different stellar masses and ages and under- 
stand the evolution of disk structures and grain 
growth, and therefore the environment of possi- 
ble ongoing planet formation. In this paper we 
present near- infrared coronagraphic imaging of RY 
Tau from the SEEDS program, the first publica- 
tion of the near-infrared scattered light associated 
with the disk around this star. RY Tau is an ac- 
tive pre-main sequence star with a stellar mass of 
2 M 0 (Calvet et al. 2004; Isella et al. 2009). The 
estimated age ranges from 0.5 (Isella et al. 2009) to 
8 Myr (Calvet et al. 2004). The star is associated 
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with a relatively massive disk (3x 10 -3 — 10“ 1 Mq, 
Isella et al. 2009, 2010) with a large infrared ex- 
cess at near-to-far infrared wavelengths (e.g., Ro- 
bitaille et al. 2007), an optical jet (St-Onge and 
Bastien 2008; Agra-Amboage et al. 2009), a scat- 
tering nebulosity due to the remnant of the en- 
velope (e.g., Nakajima and Golimowski 1995; St- 
Onge and Bastien 2008), and a large time varia- 
tion in optical photometry and spectroscopy (e.g., 
Petrov et al. 1999). These indicate a relatively 
young evolutionary phase. Despite this, the H<a 
equivalent width observed over the last 30 years is 
relatively low (8-20 A, see Chou et al. 2013, and 
references therein), similar to more evolved pre- 
main sequence stars. 

Recent observations with millimeter interferom- 
etry seem to show evidence of a hole in the disk 
seen at 1.3-mm, with a radius of 15 AU (Isella 
et al. 2010). Such a hole can be made by tidal in- 
teraction between the inner disk and protoplanets 
(e.g., Papaloizou 2007; Zhu et al. 2011). These 
disks are called transitional disks (e.g., Hughes 
et al. 2009; Brown et al. 2009, and references 
therein; Mayama et al. 2012). However, RY Tau 
is different from the other transitional disks be- 
cause it still drives a jet and has no evidence for a 
deficit of warm thermal emission near 10 /im (e.g., 
Robitaille et al. 2007). 

The rest of the paper is organized as follows. 
In Section 2 we summarize our observations and 
data reduction. In Section 3 we show the ob- 
served PI distribution and polarization. We then 
perform comparisons with simulations of scattered 
light using conventional disk and dust models with 
the two approaches described below. In Section 
4 we use the full radiative transfer tools devel- 
oped by Whitney et al. (2003b); Robitaille et al. 
(2006, 2007) with a disk geometry obtained by fit- 
ting the spectral energy distribution (SED). The 
model SEDs include the processes of scattering, 
absorption, and re-radiation on dust grains at all 
the wavelengths from UV to radio. In Section 5 we 
use monochromatic simulations of scattered light 
with our own dedicated code to attempt to bet- 
ter fit the observed PI distribution. Although 
the thermal structure and re-radiation process are 
not included, the emission from the disk surface is 
dominated by scattered light on dust grains, and 
this simplification allows us to conduct simulations 
covering a large parameter space in the disk geom- 


etry and grain size distribution. 

In Section 6 we discuss the implications for the 
scattering geometry and dust grains, and the pos- 
sible origins of the non-axisymmetry in the ob- 
served PI distribution. Throughout the paper we 
adopt a distance to the target of 140 pc (Wich- 
mann et al. 1998). 

2. Observations and Data Reduction 

Observations were made on 2011 January 27 
using HiCIAO and AO 188 at Subaru 8.2-m. As 
with several other SEEDS observations, the polar- 
ization was measured by rotating the half wave- 
plate to four angular positions (in the order 0°, 
45°, 22.5°, and 67.5°) using the PDI (Polariza- 
tion Differential Imaging) +ADI (Angular Differ- 
ential Imaging) mode. A single Wollaston prism 
was used to split incident light into two images, 
each with a 20” x 9” field of view and a pixel scale 
of 9.5 mas pixel -1 . We obtained 13 full waveplate 
rotation cycles, taking a 30-s exposure per wave- 
plate position, with a 0”. 3-diameter coronagraphic 
mask. The field rotation was ^8° during the ob- 
servations. 

The data were reduced using the standard ap- 
proach for polar imetric differential imaging (Hink- 
ley et al. 2009) as well as the other SEEDS studies 
(e.g., Hashimoto et al. 2011, 2012; Kusakabe et al. 
2012). The reduction was made using the Image 
Reduction and Analysis Facility (IRAF) 2 , pyRAF 
and python. Telescope and instrument polariza- 
tion was corrected following Joos et al. (2008). 

We find that the intensity (/) distribution ob- 
served in the disk region varies between expo- 
sures. This implies that the correction of the 
point-spread function with adaptive optics (AO) 
was not stable during the observations due to the 
moderate quality observing conditions, resulting 
in the stellar flux leaking out to a different de- 
gree in different images. We therefore use four im- 
ages (at the waveplate positions of 0°, 45°, 22.5°, 
and 67.5°) with the minimum intensity distribu- 
tion (i.e., minimum flux for the halo of the point- 
spread function associated with the star) to derive 
a lower limit for the degree of polarization. Even 


2 IRAF is distributed by National Optical Astronomy Obser- 
vatory, which is operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agree- 
ment with the National Science Foundation. 
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in these images, the / distribution is centrosym- 
metric, i.e., very different from that expected for 
the disk associated with RY Tau (Sections 3-5), 
suggesting that they are severely contaminated by 
the stellar flux. 

The unstable AO correction also resulted in dif- 
ferent PI distributions during the observations. 
To investigate this effect, we have calculated the 
co-added PI images in a few different ways: i.e., 
taking an average or median for all the data set 
(i.e., 13 full waveplate rotation cycles), or select- 
ing the best data sets and averaging them. All of 
the methods provide almost identical results. In 
the rest of the paper, we use the median images 
for the 13 full waveplate rotation cycles, with a 
total integration time of 1560 s. 

Before obtaining the coronagraphic frames we 
observed the object without a coronagraphic mask 
to measure the integrated stellar I flux (hereafter 
/*) and normalize the PI flux for each pixel. We 
obtained three frames with a 1.5-s exposure, with 
a half-wave plate P.A. of 0° and an l-% ND filter. 
The flux /* is obtained by integrating the o- and e- 
fluxes over the space in the same exposure. From 
these images and the individual science images we 
also estimate a median Strehl ratio of 0.30-0.37 
during the RY Tau observations. 

After observing RY Tau with the coronagraphic 
mask we observed a reference main sequence star 
(HD 282411) with three 10-s exposures, with a 
half-wave plate P.A. of 0° and a 0.1-% ND filter, 
and without a coronagraphic mask. These images 
(average Strehl ratio of 0.56) will be used to con- 
volve the simulated images in the later sections. 

From the above images we derive an H- 
magnitude for RY Tau of 5.7. This is ^0.4 mag 
brighter than the previous measurements of 6.1 
mag by Kenyon and Hartmann (1995) and the 
2MASS all-sky survey (e.g., Robitaille et al. 
2007). This discrepancy can be attributed to vari- 
ability due to obscuration by circumstellar dust 
(Sections 6.3). 

3. Observed PI Distribution 

Figure 1 shows the distribution of the PI flux 
overlaid with the polarization vectors. The PI 
distribution is elongated along the major axis of 
the disk observed via millimeter interferometry 
(P.A.=24°/204°, Isella et al. 2010). Its angular 


scale is ^1”.0 (^140 AU) and ^0”.6 (^80 AU) 
along the major and minor axes, respectively. The 
angular scale for the major axis is similar to that 
observed by Isella et al. (2010). 

The bright part of the PI emission is offset from 
the star towards the blueshifted jet (P.A. =294°, 
St-Onge and Bastien 2008; Agra-Amboage et al. 
2009), i.e., the far side of the disk. The PI distri- 
bution shows a minimum along the direction of the 
jet (P.A. =294°), and increases to maxima at P.A.s 
of ^210° and ^350°. This butterfly-like morphol- 
ogy in the PI distribution is similar to that mod- 
eled for some of the flared disks seen at optical 
wavelengths by Min et al. (2012). While the mor- 
phology in the PI distribution is relatively sym- 
metric about the jet (and disk) axis, its brightness 
is asymmetric. The southwest side is brighter than 
the northern side by a factor of ^2 near the peaks 
(PI/P ~ 8 x 10 -7 and ~ 3 x 10 -7 per pixel, re- 
spectively). 

The observed PI distribution, asymmetric with 
respect to the major axis of the disk, contrasts 
with the other disks observed in the SEEDS pro- 
gram (Hashimoto et al. 2011, 2012; Muto et al. 
2012; Kusakabe et al. 2012; Tanii et al. 2012; 
Mayama et al. 2012; Grady et al. 2013; Follette 
et al. 2013), in which the PI distribution is gener- 
ally symmetric. This indicates that RY Tau is as- 
sociated with a geometrically thick disk or a rem- 
nant envelope, and therefore at a younger evolu- 
tionary stage than the others (see, e.g., Fischer 
et al. 1994; Whitney et al. 2003a). This agrees 
with the stellar age of ^0.5 Myr estimated by 
Isella et al. (2010), i.e., the star is in the earliest 
stage of the Class II phase, in which the star be- 
comes visible at optical wavelengths (e.g., Stahler 
and Palla 2005). 

Simulations show that such a disk is associated 
with faint emission at the other side of the disk 
(Fisher et al. 1994; Whitney et al. 2003a; Sec- 
tion 5). Although Figure 1 shows a similar PI 
distribution in the northeast to south of the star, 
this region is contaminated by an artifact caused 
by the variation of the halo in the point-spread 
function (PSF) of the star during the modest ob- 
serving conditions. From our image we estimate 
an upper limit for PI/P of 0.4 x 10 -7 per pixel. 

The polarization vectors show a centrosymmet- 
ric pattern, as observed in several other disks. 
This is consistent with the scattered flux being 
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dominated by a single scattering, with multiple 
scatterings being negligible (Sections 6.2). As de- 
scribed in Section 2, the degree of polarization 
shown here is a lower limit due to contamination 
from the PSF halo of the star in the / flux distri- 
bution. We measure a lower limit for the degree 
of polarization of 9 db 1 % at the peak in the south- 
west. The actual degree of polarization due to the 
scattered light would be significantly higher than 
this limit. 

4. Conventional Full- Radiative Transfer 
Models with Fitting Spectral Energy 
Distributions 

Robitaille et al. (2007) developed a fitting tool 
for the SEDs of young stellar objects, and applied 
it to sources including RY Tau to investigate the 
physical parameters of the star, disk and enve- 
lope. This fitting is made via comparison between 
the observed SED and a grid of 200,000 modeled 
SEDs covering a wide range of parameter space 
calculated by Robitaille et al. (2006). Once we 
derive the parameters for the star, disk, and en- 
velope using the above SED fitter, we can obtain 
the simulated PI image of the scattered light in 
the disk and envelope using the Robitaille et al. 
(2006) code. 

In this section we use this approach to attempt 
to reproduce the observed PI image. In Section 
4.1 we briefly summarize the models based on the 
description in Robitaille et al. (2006); Whitney 
et al. (2003a, b), and our procedure for obtaining 
the simulated PI images. In Section 4.2 we show 
the results of the SED fitting and the PI image of 
the disks based on the fitted parameters. 

Note that we use an axisymmetric distribution 
of the circumstellar material for the simulations in 
this and following sections. This implies that, in 
principle, these simulations cannot reproduce the 
asymmetry at the PI distribution about the rota- 
tion axis of the disk (and the jet axis) described 
in Section 3. Simulations for an asymmetric PI 
distribution are beyond the scope of this paper. 
The implications for this asymmetry is discussed 
in Section 6.3. 

4.1. Robitaille et al. (2006, 2007) Models 

The modeled system consists of a pre-main se- 
quence star with an axisymmetric circumstellar 


disk, an infalling flattened envelope, and an out- 
flow cavity. The SED is determined using full ra- 
diative transfer in the disk and envelope, i.e., in- 
cluding absorption, scattering and re-radiation of 
light by dust grains in addition to the flux directly 
observed from the star. The heating sources of the 
disks are stellar radiation and viscosity in the disk. 
The stellar spectrum is parameterized by the star’s 
effective temperature and radius, via comparisons 
with modeled spectra for stellar photospheres (Ku- 
rucz 1993; Brott and Hauschildt 2005). 

The density distribution of a standard flared 
accretion disk (e.g., Shakura and Sunyaev 1973; 
Lynden-Bell and Pringle 1974) is described in 
cylindrical coordinate (r,z) by: 


P(r , z) = po 1 



.h. 


(i) 

where po is a constant to scale the density; R* 
is the stellar radius; a is the radial density expo- 
nent; and h is the disk scale height. The scale 
height h increases with radius as h = hor^, where 
/ 3 is the flaring power (/3 > 0). According to 
more detailed models of disk structures by Chiang 
and Goldreich (1997); D’Alessio et al. (1999a), the 
scale height exponent varies with radius in the in- 
ner disk but follows a similar power law beyond 
a few AU (Cotera et al. 2001). In addition to 
the above parameters the minimum and maximum 
radii of the disk are included as free parameters. 
The gas-to-dust mass ratio is assumed to be 100. 


Throughout the simulations Robitaille et al. 
(2006, 2007) assume a = /3 + 1. This yields the 
surface density distribution E(r) oc r -1 , approxi- 
mately agreeing with that inferred from millime- 
ter interferometry for disks associated with many 
low-mass pre-main sequence stars (Andrews et al. 
2009, 2010). This surface density power index may 
not be consistent with the one for the RY Tau 
disk seen at millimeter wavelengths by Isella et al. 
(2010). Even so, models with the above assump- 
tion fit the observed SEDs well, as shown in Sec- 
tion 4.2. 


The density structure for the envelope is given 
by Ulrich (1976); Terebey et al. (1984). The enve- 
lope is associated with a cavity whose shape varies 
as 2 oc r 1-5 . Robitaille et al. (2006) contains fur- 
ther details for the envelope and outflow cavity. 
Although the envelope mass, cavity density etc. 
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will be calculated using the SED fitter, their con- 
tribution to the SEDs of optically visible pre-main 
sequence stars such as RY Tau will be negligible 
(Section 4.2). 

The dust composition assumes a mixture of as- 
tronomical silicates and graphite in solar abun- 
dance without an ice coating. The following two 
grain size distributions are used based on Wood 
et al. (2002); Whitney et al. (2003a, b): (1) the dis- 
tribution in the denser regions of the disk (m# 2 > 
10 10 cm -3 ) where one would expect significant 
grain growth up to 1 mm; and (2) the distribution 
in the more diffuse regions (m# 2 < 10 10 cm -3 ) in 
which the grain sizes are slightly larger than those 
in the diffuse interstellar matter. The grain size 
distribution for the former is described in Wood 
et al. (2002), and that for the latter is similar to 
Kim et al. (1994) (hereafter KMH; Whitney et al. 
2003a). The former affects the thermal balance 
in the disk via radiative transfer, and is success- 
ful in fitting the SED of the HH 30 disk (Wood 
et al. 2002; Whitney et al. 2003a, b). The latter 
grains are responsible for scattered light on the 
disk surface (and in the envelope) at optical and 
near-infrared wavelengths. The optical constants 
for silicate and graphite are taken from Laor and 
Draine (1993). 

The 200,000 SEDs provided by the Robitaille 
et al. (2006) code include flux and polarization 
spectra for 250 wavelengths (from 0.01 to 5000 
/i m). The authors used 20,000 parameter sets and 
computed the results for 10 viewing angles from 
face-on to edge-on at equal intervals in the cosine 
of the inclination. The SED fitter developed by 
Robitaille et al. (2007) searches for best fitting 
SEDs using the minimum y 2 method to fit the 
observed fluxes at a range of wavelengths. We set 
a distance to the object of 140 pc, and an accept- 
able range for the visual extinction Ay of 2. 0-2. 4, 
based on Calvet et al. (2004). See Robitaille et al. 
(2007) for further details of the fitting process. 

To derive the star /disk/envelope parameters for 
RY Tau we used their SED fitter using the pho- 
tometric data tabulated in Robitaille et al. (2007) 
(Table 1). These parameters are used to model the 
PI images with the Robitaille et al. (2006) radia- 
tive transfer code. We used 10 7 photons for each 
case. The resultant images for the Stokes parame- 
ters Q and U were convolved with the PSF of the 
reference star (Section 2) before obtaining the PI 


image using PI = \] Q 2 + U 2 . Then the PI flux is 
normalized to the stellar I flux /*, and scaled to 
match the pixel size of Subaru-HiCIAO. 

4.2. Results 

Figure 2 shows the observed and ten best-fit 
SEDs. We note that the fluxes at different wave- 
lengths were measured at different epochs (Table 
1), and may be highly time- variable. In particular, 
Herbst et al. (1994) reported a variation in the op- 
tical fluxes of a factor of ^8 from 1961-1980 (see 
also Petrov et al. 1999, for longer monitoring ob- 
servations). Even so, the fluxes in Figure 2 are well 
fitted by a single SED except for the 12 /im flux ob- 
served using the InfRared Astronomical Satellite 
(IRAS). A larger excess at this wavelength may be 
attributed to bright silicate emission (Honda et al. 
2006) in the filter coverage (8.5-15 /im). 

The parameters for the modeled SEDs are 
shown in Table 2 along with the y 2 value and 
the model ID specified by Robitaille et al. (2006). 
The constant used to scale the scale height rela- 
tion (ho) is fixed at r=50 AU from the star. The 
modeled SEDs are based on four sets of physical 
parameters for the star /disk/envelope with incli- 
nation angles of 57°-76°. The inclination angle 
for most of the SEDs and the disk mass for all the 
models are approximately consistent with those 
of Isella et al. (2010) obtained using millimeter 
interferometry (> 65° and >3x 1O _3 M 0 , respec- 
tively) . 

Despite the similarity in the shape of SEDs, 
some parameters are significantly different be- 
tween models. The outer disk radius is ^80 AU for 
two out of four models, but ^120 and ~400 AU 
for the remaining. The inner disk radius calcu- 
lated with regard to the sublimation radius varies 
by a factor of ^8 between models; the disk mass 
by a factor of ^4; disk accretion rate by a fac- 
tor of ~20; and the envelope mass by a factor of 
^ 10 4 . The gas+dust density in the outflow cav- 
ity is ^ 2 x 10 -21 g cm -3 in two out of the four 
models, but zero for the remaining. 

Figure 3 shows the simulated PI images for the 
four physical parameter sets, with different incli- 
nation angles. Despite significant differences in 
the parameters described above, the four param- 
eter sets result in PI distributions that are strik- 
ingly similar to each other, except for the fact that 
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model 3007615 produces a slightly fainter distribu- 
tion than the others, due to the small flaring of the 
disk (ft-50AU = 1-6 AU, f3 = 1.12). This similarity 
between models, both for the SEDs and PI dis- 
tribution, can be attributed to the fact that both 
the infrared radiation and the scattered light at 
the disk surface are determined by the surface ge- 
ometry of the disk. The PI distribution shown in 
Figure 3 is dominated by scattering on the disk 
surface, and the contribution from the envelope 
and outflow cavity is significantly smaller. To in- 
vestigate this, we removed the envelope from the 
models with a relatively massive envelope and dif- 
fuse dust layer in the cavity (models 3000949 and 
3012376 according to the SED fitter), and ran the 
simulations again. We find that the PI distribu- 
tion is almost identical to the model with the en- 
velope and cavity, in particular within ^100 AU 
of the star. 

The models at inclination angles i— 57° and 76° 
show two separated bright regions similar to the 
observations. Furthermore, the modeled PI flux 
normalized to the stellar I flux (PI/P) is of the 
same order as the HiCIAO observations. However, 
the offset from the major axis of the disk is signif- 
icantly smaller for ah of the modeled PI images 
compared to the observations. These offsets for 
the models and observations are clearly shown in 
the bottom right plots of Figure 3. While that of 
the observed image is clearly offset from the ma- 
jor axis, the modeled images are more symmetric 
about the major axis of the disk. This indicates 
that the disks inferred from the SEDs are geomet- 
rically thinner than that producing the observa- 
tions 3 * * , if the shape of the emission is attributed 
to the disk. 

5. Monochromatic Simulations with Con- 
ventional Disk Models 

In order to explore the effect of geometry, we 
relax the density distribution of the disk con- 
strained by the SED fitting and focus on repro- 
ducing the PI distribution. To do this we carry 
out monochromatic radiative transfer calculations 
using Monte- Carlo code developed by us specih- 


3 The “thickness” here does not imply the scale height in 

Equation (1) but that between the midplane and disk sur- 

face determined with r ~ 1 from the star (e.g., Watson 

et al. 2007, see also Section 5.1 and Figure 6) 


cally for use with SEEDS observations (the Sprout 
code). We perform simulations over a large region 
of parameter space for the disk geometry, and also 
investigate the effects of different grain size distri- 
butions. We describe the details of the simulations 
in Section 5.1, and the results in Section 5.2. 

In this section we focus on comparisons between 
the models and observations for the bright part of 
the PI distribution. Although the modeled PI 
distribution is associated with faint emission at 
the other side of the disk, we will not discuss this 
component. For this component one could easily 
reproduce a PI flux consistent with observations 
(i.e., the upper limit), by adding extinction by a 
flattened envelope not included in the models in 
this section. 

5.1. Models and Simulations 

We follow the method described in Fischer et al. 
(1994). We place a central unpolarized light source 
equivalent to the star as the starting point for cal- 
culating the scattering of photons from the disk. 
The light path for the next scattering position 
is calculated for an opacity distribution based on 
the disk in Equation (1) and the dust opacity de- 
scribed below. The scattering angle and Stokes 
parameters after scattering are calculated based 
on Mie theory. The Stokes parameters for each 
photon are initially set to (/o,Qcb^o,Vb)=(l, 0,0,0) 
and normalized to 7 out = albedo • I[ n after each 
scattering. 

We use 10 6 photons for each simulation. The 
photons escaping from the disk are collected in 
imaginary detectors at different viewing angles. 
For our purpose unscattered stellar photons are 
not collected. In order to normalize the PI flux 
to the stellar I flux, we separately calculate the 
expected number of photons for each viewing an- 
gle based on the incident number of photons and 
extinction. 

For dust grains we use homogeneous spheri- 
cal particles as commonly used in other studies, 
including Whitney and Hartmann (1992); Whit- 
ney et al. (2003a, b); Cotera et al. (2001); Wood 
et al. (2002); Robitaille et al. (2006). We use 
the grain compositions of Robitaille et al. (2006, 
2007); Cotera et al. (2001); Wood et al. (2002), 
i.e., astronomical silicate and carbon dust with- 
out an ice coating. For the size distribution, we 
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use (1) that of interstellar dust measured by KMH 
(. Rv = 3.1); (2) the larger size distributions used 
by Cotera et al. (2001) and Wood et al. (2002) to 
reproduce the scattered light observed in the HH 
30 disk (hereafter C01); and (3) same as C01 but 
with the grain size scaled by a factor of 15, preserv- 
ing the total dust mass (hereafter C01xl5) (Fig- 
ure 4). Note that the C01xl5 distribution would 
be too large for the disk surface of RY Tau. Honda 
et al. (2006) conducted mid-infrared spectroscopy 
and analysis of silicate emission toward a num- 
ber of pre-main sequence stars. These authors fit 
the spectrum of individual stars using models with 
two dust sizes (0.1 and 1.5 /i m), and showed that 
the disk surface of RY Tau has one of the small- 
est ratios of large-size grains to small-size grains of 
their sample (34 %). This suggests that the grain 
size distribution in RY Tau is fairly close to that 
of the HH 30 disk, which has a central star that 
is significantly younger than most optically visible 
pre-main sequence stars (e.g., Burrows et al. 1996; 
Watson and Stapelfeldt 2007). However, we still 
show the results obtained with C01xl5 to inves- 
tigate how the results are affected by the assumed 
size distribution. 

Different authors use different types of carbon 
dust, either graphite or amorphous carbon (Cotera 
et al. 2001; Wood et al. 2002). While graphite 
has been extensively used (e.g., Draine and Lee 
1984; Laor and Draine 1993; Kim et al. 1994; 
Whitney et al. 2003b; Robitaille et al. 2006; Dong 
et al. 2012a, b), far-infrared SEDs of young stellar 
objects and evolved stars suggest the absence of 
graphite and presence of amorphous carbon in cir- 
cumstellar dust (Jager et al. 1998, and references 
therein). We still use graphite for the KHM dis- 
tribution for consistency with the authors as their 
size distribution is determined assuming graphite 
for the carbon dust. We use amorphous carbon for 
C01 and C01 xl5 for the above reason and follow- 
ing Cotera et al. (2001); Wood et al. (2002). The 
use of different carbon dusts does not significantly 
affect the modeled PI distribution, and it does not 
affect the conclusions of the paper (Appendix A). 

Calculations for Mie scattering are made using 
the code developed by Wiscombe (1996). The 
optical constants for astronomical silicate and 
graphite are obtained from Draine and Lee (1984): 
they are identical in the near infrared with the 
Laor and Draine (1993) values used for the simu- 


lations in Section 4. For amorphous carbon we use 
the optical constants of Jager et al. (1998) with a 
pyrolysis temperature of 600° C. As for the use of 
different carbon materials (amorphous carbon or 
graphite) , amorphous carbon with different pyrol- 
ysis temperatures does not significantly affect the 
modeled PI distribution and conclusions of the 
paper (Appendix A). 

Table 3 shows the physical properties (opac- 
ity, albedo, forward throwing parameter and the 
maximum degree of polarization) for the above 
three dust models. Figure 5 shows the scatter- 
ing properties of these dust models at 1.65 /im, 
i.e., /, PI, and the degree of polarization as a 
function of scattering angle, in the case where the 
incident light is not polarized. The I and PI 
fluxes for each scattering angle are derived by nor- 
malizing the scattering matrix elements S n and 
—S 12 by a constant Iq to match the results to the 
“weighted photon method” described above (i.e., 
f I /IodQ = albedo). The forward scattering is 
more significant for large grains models, but dif- 
ferent dust models show a relatively similar dis- 
tribution of scattered intensity I/Iq at scattering 
angles of 40-180°. The polarization shows a sim- 
ilar dependency on scattering angle, peaking at 
~ 90°, but the absolute value for the C01 and 
C01xl5 models are 30-40 % lower than for the 
KMH. The PI flux shows a maximum at scatter- 
ing angles of 0= 75°, 60°, and 28° for the KMH, 
C01, and C01xl5 models, respectively, and de- 
creases towards both sides for all these dust mod- 
els. 

We assume an outer radius for the disk of 80 
AU, based on our HiCIAO observations and mil- 
limeter interferometry by Isella et al. (2010). The 
inner radius is set to 1 AU. We find that our re- 
sults are not significantly affected by the choice 
of inner radius as long as it is within the coro- 
nagraphic mask (r ^30 AU at a distance to the 
target of 140 pc) and the disk surface has a flared 
geometry. Note that the millimeter interferometry 
of Isella et al. (2010) suggests the presence of a disk 
hole with a radius of ~15 AU, while the presence 
of near-infrared emission suggests the presence of 
warm dust within 1 AU (Akeson et al. 2005; Pott 
et al. 2010). 

The remaining free parameters for the density 
distribution of the disk are po? ^soau, /?, and a. 
The first three parameters all have a large effect 
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on the offset of the PI distribution from the major 
axis of the disk. This makes it difficult to search 
for the parameter sets which best fit the results. 
To overcome this problem, we set the total optical 
thickness between the star and the edge of the 
disk to 1 at a direction 0 from the midplane of 
the disk. This angle will be adjusted to fit the 
degree of the offset of the PI distribution from 
the major axis of the disk (see Figure 6 for an 
example). This determines the parameter po> and 
as a result, our parameter searches will be made 
with the remaining two free parameters, i.e., ft, 5 oau 
and (3. For the simulation in Section 5.2 we adopt 
a = (3 + 1 as we did in Section 4. The use of an 
independent a produces almost identical results. 

In addition to po? ^soau? /?, and <a, Equation 
(1) also includes the stellar radius R*, but its un- 
certainty does not affect the results of simulations. 
The stellar radius is significantly smaller than the 
radii for the density distribution we use (i.e., >1 
AU), and therefore the term y/R*/r produces only 
a very minor contribution to Equation (1). The 
equation also includes the term (R*/r) a , but the 
stellar radius here can be replaced by any radius, 
depending on the radius at which we define po- 
Again, this parameter is not explicitly provided 
but scaled using the constraint of the optical thick- 
ness described above. 

Throughout the simulations we assume that re- 
radiation from the disk is negligible, or that it oc- 
curs in very close proximity to the star so that the 
effect of the extended structure in the emission 
region far from the star is negligible. To prove 
the validity of this assumption, we ran this simu- 
lation code with the disk parameters in Table 2, 
and found that the code reproduces the PI distri- 
bution in Figure 3 well. 

5.2. Results 

Figure 7 shows examples for the PI distribution 
using the C01 grain size distribution (i.e., that of 
the HH 30 disk) . The combination of the angle for 
r = 1 (0 = 30°) and the viewing angle (49° from 
the face-on view) are selected to approximately fit 
the observed morphology. As for the observations 
and models in Sections 3 and 4, the PI flux dis- 
played is normalized to the stellar I flux. The 
image convolution is made in the same manner as 
Section 4. The contrast for each figure is adjusted 
in each image for the best morphological compar- 


ison for the bright part. 

In Figure 7 all except the case of ft 50 AU= 25 AU, 
/3= 1.3 show a butterfly- like morphology similar to 
the observed PI distribution. The peak PI flux in 
these models, however, is 3-6 times larger than the 
observed value. Increasing f3 produces a brighter 
PI distribution. In the case of (3= 2.0 and 2.7, the 
modeled PI /I* is larger for ft50AU=25 AU than 
^50AU=5 AU by a factor of ^ 1.5. This difference 
is attributed to the fact that, for the latter model, 
the star suffers extinction from the edge of the disk 
which decreases /*. In the case of (3= 2.0 and 2.7, 
increasing ftsoAU also produces a larger PI flux at 
the other side of the disk. The PI distribution 
is fairly centrosymmetric for ft 50 AU= 25 AU and 
/3=1.3 for any contrast. 

In the same figure we also show the density pro- 
files for the individual disks. The dashed, solid, 
and dotted curves show the positions at which 
t= 0.5, 1, and 2, respectively, measured from the 
star. These curves indicate that the scattering 
layer is geometrically thin at the disk surface. In- 
creasing j3 produces a large flaring angle in the 
outer disk, which explains the brighter PI flux as 
the surface is more easily illuminated. The cur- 
vature of the surface remains similar for the same 
/?, however, the disk is not flared for /3=1.3 and 
^50 au- la this case, the scattering is dominated 
by the inner part of the disk, producing a cen- 
trosymmetric PI distribution as described above. 

The viewing angle selected for Figure 7 (49°) 
is smaller than that determined from millimeter 
emission by Isella et al. (2010) (> 65°). A larger 
viewing angle does not reproduce the observed PI 
morphology. Figure 8 shows examples for a view- 
ing angle of 63°. The disk parameters are identi- 
cal to those for Figure 7, but the optical thickness 
along the radial direction has been set to r = 1 at 
0 = 25° to fit the offset of the bright PI emission 
from the major disk axis. The results are similar 
to those of Figure 7, but the bright part shows a 
thinner distribution in the vertical direction than 
the observations. Furthermore, a large (3 results 
in a larger spatial extension at the bright side. To 
clearly show the former discrepancy, we extract 
the 1-D profiles at the positions indicated in Fig- 
ures 7 and 8 and show them in Figure 9. The peak 
PI fluxes are comparable to or larger than those 
in Figure 7, and therefore significantly larger than 
the observations. 
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Figure 10 shows the PI images for the same 
disk parameters as Figures 7 and 8 but with differ- 
ent dust models. Changing the dust model results 
in little variation in the PI images. This can be 
explained by the fact that the PI flux at the bright 
side is dominated by a single scattering with scat- 
tering angles 60° — 135°, and the PI flux for this 
range of angles is similar for the three different 
dust models, as shown in Figure 5. In particular, 
the PI fluxes shown in Figure 5 for these scatter- 
ing angles explain the fact that the PI flux for the 
KMH model is larger by a factor of ^2 than the 
others in Figure 10. Table 4 shows the peak PI/P 
for individual models. These are larger than the 
observations by a factor of 2-12 at the peak. 

In summary, the modeled morphology in the PI 
distribution does not match the observations with 
the viewing angle inferred from the millimeter ob- 
servations (>65°). As shown in Figures 7, 8, and 
10, the modeled PI distribution from the bright 
side of the disk does not vary significantly with 
the value of h$ oau> P, or the choice of dust model 
once we determine the outer radius of the disk and 
the optical thickness in a given direction. Some 
models may also be excluded due to the degree of 
polarization, the inferred disk mass or extinction 
toward the star. However, these do not provide 
constraints as clearly as those from the PI flux. 
We discuss these constraints in Appendix B. 

6. Discussion 

In Section 6.1 we discuss the possible pres- 
ence of an optically thin scattering layer above 
the disk which could be responsible for the ob- 
served PI distribution. In Section 6.2 we briefly 
discuss whether the use of more realistic dust mod- 
els might still fit the observations using a conven- 
tional disk model. In Section 6.3 we discuss the 
implications of the asymmetry of the PI distri- 
bution about the jet/disk axis that was shown in 
Section 3. 

6.1. Geometry of the Scattering Layer 

In Sections 4 and 5, we found inconsistencies 
between the models obtained with conventional 
disk and dust models and the observed PI dis- 
tribution. The disk geometry obtained from the 
SED corresponds to a PI distribution with a sig- 
nificantly smaller offset from the major axis of the 


disk (Section 4). On the other hand, the disk 
geometries which reproduce this offset require a 
smaller viewing angle (i.e. , close to face-on) than 
that determined through millimeter observations 
(Section 5). We have not explored the full pa- 
rameter space, but once we fix the parameters to 
reproduce the angular scale along the major axis 
and the offset of the bright PI emission on the far 
side of the disk, varying disk and dust parameters 
has little effect on the resulting PI distribution. 

The results of the two sets of simulations (full 
radiative transfer and monochromatic scattering) 
may imply that the system consists of (1) a ge- 
ometrically thin disk which is partially responsi- 
ble for the infrared SED but does not contribute 
to the PI flux in the near infrared; and (2) an 
optically thin and geometrically thick upper layer 
which is responsible for the observed PI distribu- 
tion in the near infrared and the remaining mid- 
to-far infrared flux. A similar geometry has been 
proposed for SEDs observed for some Herbig AeBe 
stars (Meeus et al. 2001, Group I). Follette et al. 
(2013) also discuss such a geometry for the SR 21 
disk in which the presence of a inner cavity is indi- 
cated by sub-millimeter observations but not seen 
in the scattered light in near-infrared. 

Figure 11 shows a schematic view for the disk, 
an optically thin layer, stellar radiation, scattered 
light and infrared re-radiation. A rough analytic 
estimate for the optical thickness, density, and 
mass of such a layer is described in Appendix C. 
The observed PI image does not show evidence 
for emission along the major axis without an off- 
set as modeled in Section 4. This can be explained 
if the thickness of the disk is significantly smaller 
than those used in Section 4, which provide a PI 
flux level comparable to the observations. 

Such a scattering geometry can easily explain 
the observed PI flux and offset from the major 
axis of the disk, by adjusting the column density 
and vertical distribution of the optically thin layer, 
respectively. Here we qualitatively demonstrate 
this using the Sprout code with an optically thin 
upper layer using Equation (1) and an optically 
thick disk with a geometrical thickness of zero. 
To focus on the scattered light outside the coron- 
agraphic mask, we set the density to be zero for 
disk radii within the mask (r < 28 AU). 

Figure 12 shows an example of simulations us- 
ing the KMH dust model, with a viewing angle 
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of 70°. The butterfly- like bright PI distribution 
approximately matches the observations. We mea- 
sure the modeled polarization degree to be 25 % 
at a position corresponding to the PI peak of the 
observations. This is also consistent with the ob- 
servations of ^>9 % measured in Section 3. The 
faint extended emission at the other side of the 
disk is brighter than the upper limit of the obser- 
vations by a factor of ~2. The modeled flux in 
this region could be adjusted, e.g., if we included 
extinction by a remnant envelope, or if we define 
an exponential cutoff for the outer radii of the op- 
tically thin layers like that used to reproduce the 
millimeter flux distribution of the disk (e.g., An- 
drews et al. 2009; Isella et al. 2010). 

The use of different dust models (C01, C01xl5) 
provides similar results to Figure 12, but the po- 
larization degree is significantly smaller: 13-14 % 
at the position corresponding to the PI peak of 
the observations. Such a polarization may be too 
close to the lower limit of the observations of ^9 
%. It is beyond the scope of the paper to identify 
the best dust models using the modeled degree of 
polarization. See Section 6.2 for more complex 
and realistic dust models than used in our study. 

Use of a geometrically thinner disk than Sec- 
tion 4 implies that the disk gets less stellar flux, 
thereby resulting in less infrared excess than that 
shown in Figure 2. This issue could be solved as 
follows. As suggested by Chiang and Goldreich 
(1997), an optically thin layer on the disk surface 
scatters or reemits directly to space about half the 
radiation, while the other half is scattered/emitted 
inward. The extra illumination from the optically 
thinner layer would warm up the disk and enhance 
the observed infrared flux. Although such a opti- 
cally thin layer is usually assumed to be geometri- 
cally thin (e.g., Chiang and Goldreich 1997; Dulle- 
mond et al. 2007), the same physical mechanism 
should work even in the case where the layer is 
geometrically thick. If the scattered light is domi- 
nated by the optically thin layer, this implies that 
the disk surface would be illuminated more signif- 
icantly from this layer than the direct stellar flux. 
Therefore, the presence of this layer may signifi- 
cantly affect the temperature structure of the disk. 

In addition, the thermal radiation from the op- 
tically thin layer may directly contribute to the 
infrared SED (e.g., Chiang and Goldreich 1997). 
However, the fraction compared to the disk flux 


would be significantly smaller than that of the 
scattered light in the near-infrared due to signif- 
icantly smaller dust opacities in the mid-to-far 
infrared (e.g., Wood et al. 2002; Whitney et al. 
2003a, b; Dong et al. 2012b). 

6.2. Dust Grains 

While homogeneous spherical particles are 
widely used for the dust in disks (Whitney and 
Hartmann 1992; Fischer et al. 1994; Whitney et al. 
2003a, b; Robitaille et al. 2006, 2007), dust parti- 
cles in practice appear to be aggregated (e.g., Hen- 
ning and Stognienko 1996; Dominik 2009), and 
may also be coated with ice (Malfait et al. 1999; 
Meeus et al. 2001; Honda et al. 2009). These facts 
would produce different optical properties. Min 
et al. (2012) conducted simulations for scattered 
light from the disk using both spherical and aggre- 
gated particles, and show that a lower PI flux is 
expected for the latter due to the extremely large 
forward scattering. This would allow us to better 
match the PI flux we simulated in Section 5 for 
the disk surface, as our values were significantly 
larger than the observations for many cases. 

However, the disks used in Section 5 would pro- 
duce significantly different SEDs than those shown 
in Figure 2. A large flaring allows the disk to re- 
ceive more stellar photons, which warms up the 
disk more efficiently, thereby producing a larger 
excess than the thinner disks we used in Section 
4. Furthermore, some disks used in Section 5 have 
a large extinction toward the star (Appendix B). 
This would result in a double peaked SED in the 
optical to near- infrared, and far infrared, respec- 
tively (e.g., D’Alessio et al. 1999b; Robitaille et al. 
2007). The former problem could be overcome 
if the dust particles had a smaller emissivity in 
mid-to-far infrared, or large albedo at UV-to-IR 
in order to impede the absorption of the stellar 
photons. To investigate further, one would need a 
detailed understanding of the optical properties of 
the aggregates over a wide range of wavelengths, 
plus simulations with full radiative transfer. 

Alternatively, observations of the PI distribu- 
tion at different wavelengths would give useful con- 
straints for the nature of the dust grains. We ex- 
pect that the PI flux is dominated by the photons 
with a single scattering, even in the case of opti- 
cally thick disks (such as those used in Sections 4 
and 5) in which multiple scattering occurs as well 
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as single scattering. The contribution from multi- 
ple scattering would be negligible for the reasons 
below. In the case of small grains (e.g., KMH), 
where photons are fairly isotropically scattered, 
the scattered photons have polarizations with a 
variety of P.A., canceling each other out. In the 
case of large grains (e.g., C01 and C01xl5), in 
which most of the photons are scattered forward, 
the polarization of individual photons is signifi- 
cantly reduced after the first scattering. 

To demonstrate this, we re-ran the simulations 
for the disks in Section 5, removing all the pho- 
tons which experience multiple scatterings. We 
find that the results are almost identical for the 
PI flux distribution, with a flux difference within 
10 % at the peak. This is in contrast with the / 
flux for which multiple scatterings are responsible 
for 20-40 % of the flux. Throughout, studies of 
the PI flux allow a more accurate investigation of 
the dust properties than studies of the I flux, in 
particular for optically thick disks. 

6.3. Implication of the Asymmetric PI 
Distribution 

As shown in Section 3, the observed PI distri- 
bution is asymmetric about the rotation axis of the 
disk: the southwest side is brighter than the north- 
ern side by a factor of ~ 2. A similar asymmetry 
in the distribution of the scattered light has been 
extensively studied for the HH 30 disk, i.e., a low- 
mass protostar at a similar evolutionary stage, us- 
ing the Hubble Space Telescope (e.g., Burrows et al. 
1996; Stapelfeldt et al. 1999; Cotera et al. 2001; 
Watson and Stapelfeldt 2007). These authors at- 
tributed such a distribution to uneven illumina- 
tion by the central object. These observations 
also show that the distribution of the scattered 
light is time variable. Watson and Stapelfeldt 
(2007) summarize possible mechanisms to explain 
the asymmetric distribution and its time variation. 
These include (1) hot (or cool) spots on the star; 

(2) shadowing by a non-axisymmetric inner disk; 

(3) obscuration by a companion star; and (4) ob- 
scuration by a disk associated with a companion 
star. The time variability in the scattered light is 
also observed in the more evolved disk system HD 
163296 (Wisniewski et al. 2008). The authors at- 
tributed this to the asymmetric shadowing of the 
inner disk. 

Of the explanations described above, shadow- 


ing by the inner disk may be the most likely for 
the asymmetric PI distribution in RY Tau. While 
uneven illumination by hot (or cool) spots is one of 
the favored explanations for the HH 30 disk (Wood 
and Whitney 1998; Stapelfeldt et al. 1999; Cotera 
et al. 2001; Watson and Stapelfeldt 2007), it is not 
likely for RY Tau: the optical continuum associ- 
ated with this star does not show clear evidence 
for time variability due to such spots (e.g., Petrov 
et al. 1999; Chou et al. 2013). Obscuration by 
a secondary star/disk may also be excluded, since 
there is no clear evidence for the presence of a close 
companion associated with RY Tau. Near-infrared 
interferometric observations have ruled out the 
presence of a companion at 0.35-4 AU from the 
star and down to a stellar flux ratio of 0.05 (Pott 
et al. 2010). 

The shadowed-disk scenario is corroborated by 
the UXOR-type variability observed toward this 
star, i.e., time- variable obscuration by circumstel- 
lar dust (e.g., Herbst et al. 1994; Petrov et al. 
1999). While many authors attribute this obscu- 
ration to the disk, its nature is not clear. It could 
be clumpy accretion of dust+gas onto the star (see 
Herbst et al. 1994; Grady et al. 2000, and refer- 
ences therein), the puffed- up inner rim of the disk 
(Dullemond et al. 2003), or the outer edge of the 
disk (The 1994). Herbst et al. (1994) also state 
the possibility that such a dust screen is associ- 
ated with a wind, not a disk or disk accretion. Of 
the above explanations, the outer edge of the disk 
may be excluded if the disk associated with RY 
Tau is geometrically thin, and we are not observ- 
ing this target close to an edge-on view (Sections 
4 and 6.1). 

RY Tau showed two abrupt brightening events 
at optical wavelengths in 1983/1984 and 1996/1997 
(e.g., Herbst et al. 1994; Petrov et al. 1999). If we 
attribute this to a single orbital period around a 
2 Mq star, the corresponding radius of the disk 
is 7 AU. This radius is comparable to the radius 
where Isella et al. (2010) identified emission peaks 
at two sides of the disk using millimeter interfer- 
ometry. If obscuration occurs in the inner disk, it 
would require time variation of the disk structures 
including warping and precession. This may be 
possible via tidal interaction with a very low-mass 
companion such as a proto-planet (Hughes et al. 
2009). 

Synoptic observations of the scattered light of 
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RY Tau would allow us to identify the origin of 
the asymmetric PI distribution, determine a typ- 
ical radius where the obscuration occurs, and con- 
strain the disk geometry (and its time variation) 
within the coronagraphic mask. Bastien (1982) 
made monitoring observations of optical polariza- 
tion integrated over the object, revealing variation 
of the position angle of polarization between -25° 
and 45° over a few years 4 . Spatially-resolved ob- 
servations like ours would have significant advan- 
tages for searching for a periodic variability toward 
this active and complicated object. We also note 
that the coronagraphic observations for the opti- 
cal I image of RY Tau were made in 2007-2009, 
and these show a brighter lobe in the north than 
the south-west (McCleary et ah, unpublished) in 
the outer envelope, i.e., the opposite side as we 
observed in the PI distribution. Although these 
observations show scattered light in significantly 
more distant outer regions, this is consistent with 
the idea of time variation of scattering light asso- 
ciated with this object. 

7. Conclusions 

We present near-infrared coronagraphic imag- 
ing polarimetry of RY Tau. The scattered light 
in the circumstellar environment was imaged at 
77-band with a high resolution (~0”.05) using 
Subaru-HiCIAO. The observed PI distribution 
shows an angular scale of ~1”.0 (~140 AU) and 
^0”.6 (^80 AU) along the major and minor axes 
of the disk, respectively, exhibiting a butterfly- like 
shape. The angular scale along the major axis is 
similar to that of the disk measured using millime- 
ter interferometry. The bright part of the emis- 
sion is offset from the star toward the direction 
of the blueshifted jet. Such a distribution can be 
explained if the object is associated with a geo- 
metrically thick disk or a remnant envelope. This 
agrees with the premise that the system is at the 
earliest stage of the Class II evolutionary phase. 

We perform comparisons between the observed 
PI distribution and simulations of scattered light 
with conventional disk and dust models. The sim- 
ulated images are made using two different ap- 
proaches. The first approach is to conduct full ra- 


4 The aperture size used for this study ranges between 8 and 
20 arcsec, and this fact may also cause different polarization 
angles. 


diative transfer simulations including absorption, 
scattering and re-radiation by dust grains, with 
disk parameters based on SED fitting (Robitaille 
et al. 2006, 2007). The second is monochromatic 
simulations with absorption and scattering only 
for a significantly larger sets of disk parameters to 
attempt to better fit the observed PI distribution. 

The first approach reproduces the PI flux level 
normalized to the integrated stellar / flux well, but 
fails to reproduce the offset of the bright PI distri- 
bution from the star along the disk axis. The sec- 
ond approach reproduces the butterfly- like mor- 
phology in PI distribution well, with a total op- 
tical thickness of 1 at 30° from the midplane of 
the disk and the viewing angle of ^ 50°. However, 
the model PI distribution does not match the ob- 
servations with a large viewing angle inferred by 
millimeter interferometry (> 65°). These results 
for the second approach are relatively independent 
of the free parameters for the disks and the size 
distribution of dust grains. 

Throughout, we find disagreements between the 
observed PI distribution and models using con- 
ventional disk and dust models. This may imply 
that the system consists of the following: (1) a 
geometrically thin disk which is partially respon- 
sible for the infrared SED but not the PI flux in 
the near infrared; and (2) an optically thin and ge- 
ometrically thick upper layer which is responsible 
for the observed PI distribution in near infrared 
and the remaining infrared flux. Simulations show 
that this idea approximately explains the observed 
PI emission with a viewing angle consistent with 
the observations. Half of the scattered light and 
thermal radiation in this layer illuminates the disk 
surface, and this process may significantly affect 
the thermal structure of the disk. 

The PI brightness has an asymmetry about the 
jet axis by a factor of ^2 in flux presumably due 
to uneven illumination caused by obscuration by 
the dusty environment. Such obscuration could 
be either due to the accretion of dust+gas onto 
the star, the puffed-up inner rim of the disk, or 
the outer edge of the disk. Synoptic studies of 
the PI distribution would give useful constraints 
for the geometry of the disk and perhaps its time 
variation within 5 AU of the star, and a better un- 
derstanding of the nature of objects with a similar 
variability at optical wavelengths (UXORs). 
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A. Simulated results with different carbon dust models 

Table 5 shows the dust properties of the C01 and C01xl5 size distributions (Section 5.1) with different 
carbon dusts: i.e. , amorphous carbon with pyrolysis temperature 400-1000 °C and graphite. The optical 
constants for the amorphous carbon and graphite are taken from Jager et al. (1998) and Draine and Lee 
(1984), respectively. Figure 13 shows their scattering properties as in Figure 5, but using the above carbon 
dusts. The simulated PI images using these dust grains show nearly identical morphologies, with differences 
in PI flux of -20 to +30 % compared with the amorphous carbon used in the main text (pyrolysis temperature 
= 600 °C). As discussed in Section 5.2, these similarities and differences are attributed to (PI/Iq) at 
scattering angles 60-135° in Figure 13. 


B. Polarization, disk masses and extinction toward the star for geometrically thick disk mod- 
els 

Figure 14 shows polarization vectors for the same modeled parameters as Figure 10. For all the models 
the vector pattern is centrosymmetric about the star as it is in the observations (Section 3). We measure a 
degree of the polarization of 10-28 % at the position corresponding to the PI peak of the observations (see 
Figure 10 for the position). Some seem inconsistent with the observations (^>9 %, Section 3). 

Tables 6 and 7 show the disk masses and extinction {Ay) toward the star, respectively, for the models 
used in Section 5.2. The disk masses are derived assuming a gas-to-dust mass ratio of 100. The extinctions 
are calculated based on the optical thickness at 1.65 /im and the ft e xt;i.65/im/ftext;0.55/im tabulated in Table 
3. Table 6 show that the disk mass is significantly smaller for H^au = 25 AU than h^AU — 5 AU with 
a given combination of the dust model, /?, and 0 T =i. As a result, the former disks produce brighter PI 
flux from the other side of the disk (Figures 7, 8, 10). Disks with a large scale height also provide a larger 
extinction toward the star. 

All the disk masses derived for H^au = 25 AU are significantly smaller than that estimated by Isella 
et al. (2010) based on interferometry (> 3 x 1O -3 M 0 ). In contrast, a majority of the disk masses derived 
for h$oAU = 6 AU are significantly larger than those of pre-main sequence stars suggested by millimeter 
observations (< 0.1 M Q , e.g., Robitaille et al. 2007; Isella et al. 2009; Williams and Cieza 2011, see also Table 
1). However, we do not constrain the models described in Section 5.2 for the following reasons: (1) the disk 
mass inferred by millimeter interferometry highly depends on the dust model used to convert the flux to 
a dust mass (Isella et al. 2010); (2) the modeled disk mass highly depends on the exponent of the radial 
density distribution <a, which is assumed to be (3 + 1 for the simulations in Section 5.2; and (3) the scattered 
light results from a small fraction to the total dust mass (Appendix C), therefore our observations do not 
directly probe the total disk mass. 

Regarding the extinction to the star, some models provide a larger value than that measured by Calvet 
et al. (2004) {Ay =2.2). However, we note that the measurement of extinction by Calvet et al. (2004) is based 
on observations at UV-optical wavelengths, adopting an extinction law for cold molecular clouds. Therefore, 
it may not be directly comparable with the extinctions of Table 7, in particular for the C01 and C01xl5 
size distributions. 


C. Optical thickness, density and mass of the optically thin and geometrically thick scattering 
layer 

Suppose the dust corresponding to the mass m is located at distance r from the star. The number of PI 
photons observed at the telescope per second is described as follows: 


Tipi — K'ext'W' 


N* 

iirr 2 



(Cl) 
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where K ext is the opacity; TV* is number of the photons at 1.65 /im ejected from the star in all directions; 
and (P/// 0 ) is the fraction of the PI flux normalized to the incident flux on dust grains (Figure 4). Q is 
the solid angle corresponding to the telescope area, thus U = A te i/d 2 , where A te i and d are the area of the 
telescope and the distance to the target, respectively. If we normalize Equation (2) by the number of stellar 
photons observed per second, i.e., n* = N*(A te i/47rd 2 ), the equation is as follows:- 


npi _ K ext m ( PA 
n* r 2 \/ 0 / ’ 


(C2) 


Here we substitute np//n* = 8 x 10 -7 , measured at the peak of the PI flux as described in Section 3; 
K ext = 1.1 x 10 4 cm 2 g -1 (C01 dust, Table 3), r=40 AU, and P/// 0 = 5 x 10 -3 str -1 (Figure 5). We 
derive a dust mass of 5 x 10 21 g. The PI flux measured at each position is based on the HiCIAO pixel 
scale (9.48 x 10 -3 arcsec), thus this dust mass corresponds to an optical thickness along the line of sight of 
r = K ex tm/A pix r^j 0.15, where A pix is the area corresponding to each pixel, assuming a distance d=140 pc. 
Therefore, it is likely that the scattering layer is optically thin along the line of sight. 

Such a layer could also be optically thin in the radial direction. The optical thickness along this direction 
highly depends on the assumed geometrical thickness of the layer. Let us assume a thickness of 30 AU, 
comparable to the thickness of the disks we used in Section 5 to explain the offset of the PI distribution 
from the star. One would expect an optical thickness of ~ 6 x 10 3 towards a HiCIAO pixel (1.3 AU) 
assuming all of the parameters described above. Integrating this over the disk radius of 80 AU, one would 
expect a total optical thickness of ^ 0.4. This is a simple estimate. In reality, the observed PI flux is lower 
than the maximum at the other positions, and the dependence of np//n* on the distance from the star 
should be included. Furthermore, the opacity at individual positions depends on the distance from the star 
(Equation 3), and this fact is not included here. 

If we assume the same layer thickness, and use the dust mass contained in the pixel with the maximum 
PI value, we would derive a dust mass density of 2.9 x 10 -20 cm -3 . This corresponds to a hydrogen number 
density of the order of ~ 10 6 cm -3 assuming a gas-to-dust mass ratio of 100. This is larger than the envelope 
density at similar radii (^50 AU) inferred by the SED fitting in Section 4 by a factor of ^10 or more. 

We estimate the dust mass of the scattering layer as follows, using Equation (3):- 

M dust = j m(r)dr = ( r PlAl (H) dr (C3) 

J J Kext \ Iq J 

where r and r are the position and the distance to the star, respectively. Substituting the projected radius 
for r, we derive a dust mass M& us t of 3 — 7 x 10 -3 M 0 . A larger value would be expected if we include 
the inclination effect for r and the smaller P///o, depending on the scattering angle at individual positions. 
Even so, this estimate would be sufficient to conclude that the optically thin layer above the disk surface has 
a mass significantly smaller than the total dust mass in the disk inferred by radio observations (10 — 50 M 0 , 
Isella et al. 2010). 
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PH I* 

(xIO 7 pixel 1 ) 



R.A. offset (arcsec) 


Fig. 1. — Distribution of PI flux and polariza- 
tion vectors in i7-band (1.65 /im). The PI flux at 
each pixel is normalized to the integrated stellar 
/ flux observed without the coronagraphic mask 
(PI /I* = 10 -7 corresponds to 6.1 mJy arcsec -2 ). 
In the PI image we set a software aperture for the 
coronagraphic mask of 0”.4 in diameter, slightly 
larger than that in the optics (0”.3 in diameter) 
to show the distribution only where the measure- 
ments are reliable. The thin dashed line shows the 
direction of the extended jet observed by Agra- 
Amboage et al. (2009) (P.A.=294°). The faint 
extended component in the northeast to south 
is probably due to an artifact due to unstable 
AO correction in the modest observing conditions. 
The degree of polarization shown here is a lower 
limit (see text). This is measured for individual 
11 x 11-pixel bins (corresponding to 0”.lx0”.l), 
and shown for those in which the mean PI /I* is 
larger than 0.4 x 10 -7 . 



Fig. 2. — Observed and modeled SEDs. Cir- 
cles are the observed fluxes tabulated in Robitaille 
et al. (2007). The error bars for the observed 
SED are shown only for those larger than the cir- 
cles. The solid curve shows the best-fit model us- 
ing their online SED fitter. The gray curves are 
the same but for the next nine best fits. See Table 
2 for the parameters for the star, disk, envelope, 
and inclination angle. 


KMH (interstellar) 

- - C01 (HH 30) 

C01x15 



a (pm) 


Fig. 4. — Size distributions of silicate and carbon 
grains for (1) interstellar medium (KMH, Kim 
et al. 1994), (2) that used for explaining the scat- 
tered light in the HH 30 disk (C01 Cotera et al. 
2001; Wood et al. 2002), and (3) that 15 times 
larger than (2) (C01xl5). 
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Model ID: 


3000949 3007615 3012376 3012808 



Disk 

inclination 

z=5 7° 


i= 76° 


Observed 



1-D Profiles ( observed 


modeled) 



100 AU 


Offset from the major axis of the disk (AU) 


Fig. 3. — Modeled PI distributions for the best-fit SED models using the full radiative transfer code de- 
veloped by Whitney et al. (2003a, b). The observed PI distribution is also shown at the bottom-left and 
in contour in the modeled images. The color contrast is the same for all the images. The contour levels 
are 0.75, 1.5, 3, and 6xl0 -7 per pixel scale of HiCIAO relative to the integrated stellar I flux, and they 
are rotated by 72° from the observations. The dashed lines in the PI images show the positions where we 
extract the 1-D distribution and show in the bottom-right. Each profile is made by averaging those at two 
sides, normalized to the peak PI flux. The observed 1-D profile is dashed for the faint side of the disk where 
it is contaminated by an artifact (Section 3). 


23 




scattering angle (deg.) 


Fig. 5. — Scattering properties of an unpolarized photon for different size distributions. From left to right 
the figures show the scattered I flux per steradian, the polarization, and the scattered PI flux per steradian, 
as a function of scattering angle. The / and PI fluxes are normalized to the incident Stokes I parameter 
(Jo), and as a result, J I /Iq dQ = scattering albedo. 
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Fig. 7. — Same as Figure 6 but for r = 1 at 6 = 30° with different h$ oau and ft. White dashed lines in the 
top-right PI image show the positions where we extract the 1-D distribution shown in Figure 9. 
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r=1 @6=25°, viewing angle=63° 
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Fig. 8. — Same as Figure 7 but for r = 1 at 9 = 25° and a viewing angle of 63° from the face-on view. 
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Fig. 10. — PI distributions with different dust models for / 3 = 2.0. The other parameters (#, viewing angle, 
and hsoAu) and contours are the same as Figures 7 and 8. 
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Fig. 13. — Same as Figure 5 but with the C01 (upper) and C01xl5 size distributions (lower) and different 
carbon dusts (amorphous carbon with pyrolysis temperatures 400/600/800/1000 °C and graphite). The 
black and gray curves are for amorphous carbon and graphite, respectively. Those for amorphous carbon 
with different pyrolysis temperatures are shown with different line styles, and that used in Section 5 (pyrolysis 
temperature of 600 °C) is shown with black thick curves. 
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r=1 @6>=30° 
viewing angle=49° 

du f t , /7 5 oau:5AU 

model 


r=1 @ 0 = 25 ° 
viewing angle=63° 



Fig. 14. — Same as Figure 14 but with the polarization vectors and without the contour for the observed 
PI distribution. 
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Fig. 6. — Density distribution of dust (left) 

and monochromatic simulations of the PI dis- 
tribution (right) for r = 1 at 0 = 20/30/40° 
(top/middle/bottom), h50AU=15 AU, (3= 2.0, 
viewing angle of 49° from the pole, and the C01 
dust size distribution. The color scale is arbitrar- 
ily adjusted to clearly show distributions similar 
to the contours, which show the observed PI dis- 
tribution. The contour scales are 0.075, 0.15, 0.3 
and 0.6 xIO -6 per pixel scale of HiCIAO relative 
to the stellar flux, and they are rotated by 72° 
from the observations. Dashed, solid and dotted 
curves show the positions for t= 0.5, 1, and 2 from 
the star, respectively. 


t=1@(9=30° t=1 @ 6 = 25 ° 

viewing angle=49° viewing angle=63° 

observed modeled 



Fig. 9. — The 1-D PI distribution extracted 
from the observed (black) and modeled PI images 
(gray). The positions of the extraction are shown 
in Figures 7 and 8. Solid and dashed gray curves 
are models for /i50AU=5 and 25 AU, respectively. 
Each profile is made by averaging those at the two 
sides and normalizing to the peak PI flux. 


NIR scattering 



Fig. 11. — A schematic view of the disk, optically 
thin layer, stellar radiation, and the near-infrared 
(NIR) scattered light and infrared re-radiation ob- 
served. See text for details. 
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Fig. 12. — A model with a disk with zero-thickness 
and an optically thin and geometrically thick up- 
per layer using the KMH dust model, (left) Den- 
sity distribution of the optically thin layer (r= 28- 
80 AU). This is set using Equation (1) with <a=2.5, 
/3= 1.5, h50AC/=30 AU, and a total optical thick- 
ness of r = 0.04 at 1° above the midplane. The to- 
tal dust mass of the layer is 7.3 x 10 3 M®. (right) 
Modeled PI distribution (color) and polarization 
vectors (lines) with a viewing angle of 70°. 3 x 10 6 
photons were used for the simulation. The color 
contrast for the PI distribution is same contrast 
as Figure 1. The black contours show the observed 
PI distribution. 
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Table 1: UV-to- Ra dio Fluxes a 


Wavelength 

(Mm) 

Flux 

(mJy) 

Year of 
Observations 

Reference 6 

0.36 

43 ± 19 b 

1962- 1990 1 ' 

1 

0.44 

(1.5 ±0.6) x 10 2b 

1962-1990 b 

1 

0.55 

(3.5 ± 1.5) x 10 2b 

1962-1990 b 

1 

0.64 

(5.1 ±2.2) x 10 2b 

1962-1990 b 

1 

0.71 

(8.9 ± 3.9) x 10 2b 

1962-1990 b 

1 

1.25 

(2.26 ±0.07) x 10 3 

1997-2001° 

2 

1.65 

(3.8 ±0.2) x 10 3 

1997-2001° 

2 

2.2 

(4.7 ±0.1) x 10 3 

1997-2001° 

2 

3.6 

(5.75 ± 0.07) x 10 3 

2003-2007 d 

3 

4.5 

(5.30 ± 0.06) x 10 3 

2003-2007 d 

3 

5.8 

(4.4 ± 0.2) x 10 3 

2003-2007 d 

3 

8.0 

(5.8 ±0.2) x 10 3 

2003-2007 d 

3 

12 

(1.774 ±0.003) x 10 4 

1983 

4 

24 

(1.8 ±0.4) x 10 4 

2003-2007 d 

3 

25 

(2.648 ± 0.005) x 10 4 

1983 

4 

60 

(1.891 ±0.007) x 10 4 

1983 

4 

70 

(9.6 ± 1.0) x 10 3 

2003-2007 d 

3 

100 

(1.4 ±0.2) x 10 4 

1983 

4 

350 

(2.4 ±0.3) x 10 3 

1990 

5 

450 

(1.9 ±0.2) x 10 3 

1990 

5 

800 

(5.6 ±0.3) x 10 2 

1990 

5 


a Cited from Robitaille et al. (2007). The L- (3.5 p, m), M- (4.8 fi m), and iV-band (10.5 p, m) data quoted from Kenyon and 
Hartmann (1995) are not included, as these have large uncertainties, and are consistent with other observations at similar 
wavelengths. 

b Measurements at multiple epochs. The uncertainty is based on the photometric variability. 

c The specific year is not clear. We describe the years of operation for the 2MASS all-sky survey. 

d The specific year is not clear. We describe the period between the launch of the Spitzer Space Telescope and when Robitaille 
et al. (2007) was published. 

e l — Herbst et al. (1994), and converted to the tabulated wavelengths by Robitaille et al. (2007) ; 2 — 2MASS all-sky survey; 3 
— Spitzer Space Telescope Archive; 4 — Weaver and Jones (1992); 5 — Andrews and Williams (2005), who quoted the values 
from Mannings and Emerson (1994) 
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Table 2 

Best-Fit Models using the SED Fitter 
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Table 3: Physical Parameters for Dust Grains (A = 1.65/im) 


Size distribution 

Mass fraction of 
carbon dust 

ftext 13 
(cm 2 g _1 ) 

^ext/^ext(0.55/i,m) 

albedo 

9 

-Fmax 

KMH (interstellar medium) 

0.38 

5.3 x 10 3 

0.16 

0.39 

0.28 

0.59 

C01 (the HH 30 disk) 

0.49 

1.1 x 10 4 

0.28 

0.45 

0.54 

0.42 

C01X15 

0.49 

5.4 x 10 3 

0.50 

0.47 

0.66 

0.37 


a See text for details. The minimum and maximum particle radii are 3 X 10 3 and 2 fj , m for KMH; 1 x 10 3 and 20 /im for C01; 
and 15 times larger for C01xl5. 

b We adopt the mass density 3.3 and 2.26 g cm -3 for silicate and graphite, respectively (Kim et al. 1994), and 1.67 g cm -3 for 
amorphous carbon (Jager et al. 1998). This parameter is often defined per gas+dust mass (e.g., Cotera et al. 2001; Wood et al. 
2002; Whitney et al. 2003a, b; Dong et al. 2012b), but in this paper we define this per dust mass to discuss the total dust mass 
in the scattering layer in Section 6.1. 
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Table 4 
Maximum PI /P 


Dust 

P 

T 

= 1 at 0 

II 

OO 

o 

o 

T 

= 1 at e = 25° 

^50AU = 

=5 AU 

25 AU 

5 AU 

25 AU 

KMH 

1.3 

3.0 X 

icr 6 

— 

3.5 x 

10 -6 

— 


2.0 

5.1 x 

icr 6 

6.5 x 10" 6 

5.0 x 

10 -6 

9.4 x 10" 6 


2.7 

4.4 x 

icr 6 

7.3 x 10 -6 

4.5 x 

10 -6 

9.3 x 10 -6 

C01 

1.3 

1.9 x 

icr 6 

— 

2.1 x 

1CT 6 

— 


2.0 

3.1 x 

icr 6 

4.0 x 10" 6 

2.7 x 

1CT 6 

5.7 x 1(T 6 


2.7 

2.6 x 

icr 2 

4.1 x 10" 6 

2.5 x 

10 -6 

6.1 x 10 -6 

C01xl5 

1.3 

1.8 x 

icr 6 

— 

1.8 x 

1CT 6 

— 


2.0 

2.6 x 

icr 6 

3.3 x KT 6 

2.3 x 

1CT 6 

5.4 x 1(T 6 


2.7 

2.1 x 

icr 6 

3.5 x KT 6 

2.0 x 

1CT 6 

5.5 x 10“ 6 


Table 5: Physical Parameters for Dust Grains with Different Carbon Dust (A = 1.65//m) 


Size distribution 

Carbon 

pyrolysis temperature 
(amorphous only, °C) 

Mass fraction of 
carbon dust 

/^ex t k 
(cm 2 g _1 ) 

albedo 

9 

-Fmax 

C01 

amorphous 

400 

0.45 

8.4 x 10 3 

0.59 

0.57 

0.37 



600 

0.49 

1.1 x 10 4 

0.44 

0.54 

0.42 



800 

0.52 

1.4 x 10 4 

0.34 

0.50 

0.48 



1000 

0.53 

1.4 x 10 4 

0.36 

0.49 

0.47 


graphite 


0.57 

1.1 x 10 4 

0.57 

0.39 

0.33 

C01xl5 

amorphous 

100 

0.45 

4.8 x 10 3 

0.58 

0.69 

0.27 



600 

0.49 

5.4 x 10 3 

0.47 

0.66 

0.37 



800 

0.52 

6.2 x 10 3 

0.41 

0.61 

0.46 



1000 

0.53 

6.0 x 10 3 

0.43 

0.60 

0.44 


graphite 

— 

0.57 

5.3 x 10 3 

0.61 

0.47 

0.24 


a See Section 5.1 and Table 3 for details. 

b We adopt the mass density 2.26 g cm -3 for graphite (Kim et al. 1994), and 1.44/1.67/1.84/1.99 g cm -3 for amorphous carbon 
with pyrolysis temperature of 400/600/800/1000 °C, respectively (Jager et al. 1998). 
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Table 6 

Disk Masses (Mq) 


Dust 

P 

r = 1 at 0 

II 

CO 

o 

o 

t = 1 at 9 = 25° 

^50AU=5 AU 

25 AU 

5 AU 

25 AU 

KMH 

1.3 

1 x 10 2 

— 

0.5 

— 


2.0 

2 

1 x HU 4 

5 x 10“ 2 

6 x 10“ 5 


2.7 

0.1 

2 x HU 4 

8 x 10“ 3 

1 x 10“ 4 

C01 

1.3 

6 x 10 1 

— 

0.3 

— 


2.0 

1 

5 x 10" 5 

2 x 10 -2 

3 x 10 -5 


2.7 

6 x 10 -2 

8 x 10“ 5 

4 x 10“ 3 

5 x 10“ 5 

C01xl5 

1.3 

1 x 10 2 

— 

5 x lO -1 

— 


2.0 

2 

1 x 10" 4 

5 x 10“ 2 

6 x 10 -5 


2.7 

0.1 

2 x 10" 4 

8 x 10 -3 

1 x 10 -4 


Table 7 

Visual extinction (A v ) toward the star 


Dust 


r = 1 at 0 

II 

CO 

o 

o 

T = 1 at 9 = 25° 

P 

^50AU=5 AU 

25 AU 

5 AU 

25 AU 

KMH 

1.3 

2 x 10“ 8 

— 

1.1 

— 


2.0 

3 x ur 7 

1.7 

1.9 

5.1 


2.7 

5 x 10" 6 

2.0 

2.6 

5.3 

C01 

1.3 

1 x 10“ 8 

— 

0.7 

— 


2.0 

2 x 10“ 7 

1.0 

1.1 

3.0 


2.7 

3 x 10“ 6 

1.2 

1.6 

3.1 

C01xl5 

1.3 

7 x 10" 9 

— 

0.4 

— 


2.0 

1 x 10“ 7 

0.5 

0.6 

1.7 


2.7 

2 x 10“ 6 

0.6 

0.9 

1.7 
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